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1 Introduction 



Hard exclusive processes are coming to the forefront of high energy nuclear and particle 
physics. This is already visible at TJNAF, COMPASS, HERMES that go for more and 
more exclusive channels, and it makes the core of the ELFE proposal. One main reason for 
that is the growing understanding that exclusive reactions provide unmatched opportunities 
to study the hadron structure, as demonstrated by recent interest to deeply virtual Compton 
scattering (DVCS) and hard diffractive meson production. All future plans also call for very 
high luminosity and would therefore be perfectly suited for the investigation of exclusive and 
semi-exclusive reactions with and without polarization. 

The classical theoretical framework for the calculation of hard exclusive processes in 
QCD was developed in 0,0], see || f|, |5j for a review. This approach introduces a concept 
of hadron distribution amplitudes as fundamental nonperturbative functions describing the 
hadron structure in rare parton configurations with a minimum number of Fock constituents 
(and at small transverse separations). The distribution amplitudes are equally important 
and to a large extent complementary to conventional parton distributions which correspond 
to one-particle probability distributions for the parton momentum fraction in an average 
configuration. The recent data from CLEO M and E791 pj provided, for the first time, 
quantitative information on the pion distribution amplitude. 

For a theorist, the main challenge is to make the QCD description of hard exclusive re- 
actions fully quantitative. Although the leading-twist QCD factorization approach correctly 
reproduces the power dependence of exclusive amplitudes on the large momentum, there 
exist many indications that soft end-point and higher-twist corrections might dominate over 
the naive QCD prediction over the range of accessible values of Q 2 . For mesons, the study of 
preasymptotic corrections to exclusive reactions has been pursued intensively and in different 
directions. In particular, meson distribution amplitudes of higher twist have been studied 
in detail in [|, ^ |9|. For baryons, the corresponding task is more complicated and received 
less attention in the past. Soft end-point contributions to several hard exclusive reactions 
involving nucleons have been estimated using various models and were found to be large [|I(J . 
To the best of our knowledge, hard higher-twist corrections and nucleon mass corrections 
have never been addressed in the literature. In this paper we present the first systematic 
study of nucleon distribution amplitudes of higher twist. 

The notion of higher twist contributions to the distribution amplitudes comprises a broad 
spectrum of effects of different physical origin: First, contributions of "bad" components in 
the wave function and in particular of components with "wrong" spin projection; second, 
contributions of transverse motion of quarks (antiquarks) in the leading twist components; 
third, contributions of higher Fock states with additional gluons and/or quark- ant iquark 
pairs. Finally, one can speak of hadron mass effects, similar to target mass corrections in deep 
inelastic scattering. The relative significance of these corrections depends on the hadron state 
in question and the particular hard process. There exists an important conceptual difference 
between mesons and baryons. For mesons, all effects due to "bad" components in the quark- 
antiquark wave functions can be rewritten in terms of higher Fock state components thanks 
to the QCD equations of motion. Since quark- ant iquark-gluon matrix elements between the 
vacuum and the meson state are numerically small (see || |§ for the existing estimates), 
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contributions of "bad" components in two-particle wave functions are small as well, and to 
a large extent dominated by meson mass corrections || ||. For baryons, on the contrary, 
QCD equations of motion are not sufficient to eliminate the higher-twist three-quark states in 
favor of the components with extra gluons, so that the former present genuine new degrees of 
freedom. Moreover, vacuum-to-baryon matrix elements of higher-twist three-quark operators 
are large, comparable to the ones of leading twist. One may expect, therefore, that higher- 
twist effects in baryon wave functions are dominated by "bad" components of three-quark 
states rather than extra gluons. A systematic study of such states as well as taking into 
account nucleon mass corrections presents the subject of this work. 

We find that a generic three-quark matrix element on the light-cone can be parametrized 
in terms of eight independent distribution amplitudes. In particular, there exists one am- 
plitude of leading twist-three that is familiar from earlier studies P, [ITfl , three distribution 
amplitudes of twist-four and twist-five, respectively, and one distribution amplitude of twist- 
six. Following the approach of || [| , we further expand all distribution amplitudes in contri- 
butions of operators with increasing conformal spin. The conformal expansion corresponds, 
physically, to the separation of longitudinal and transverse degrees of freedom, and the co- 
efficients of this expansion present the relevant nonperturbative parameters. We find that 
the higher-twist distribution amplitudes with the lowest conformal spin (asymptotic wave 
functions) involve two nonperturbative matrix elements that are well known from the QCD 



sum rule studies of the nucleon |12|, [13], [LJ]. The leading corrections to the asymptotic wave 
functions involve three more parameters (for all twists) that we estimate in this work. 

The outline of the paper is as follows. The general classification of twist-three distribution 
amplitudes is worked out in Section 2. One starts with 24 invariant functions in a general 
decomposition of a three quark light-cone operator. By isospin symmetry the number of 
independent amplitudes is reduced to eight amplitudes that determine the three-quark proton 
state completely. A convenient representation of these independent amplitudes is derived in 
terms of chiral fields. 

Section 3 is devoted to the conformal expansion. We identify the nonperturbative pa- 
rameters that describe the higher-twist distribution amplitudes up to next-to-leading order 
in the conformal expansion and work out the relations between them, as given by isospin 
constraints and equations of motion. 

In Section 4 we present particular models for the higher-twist distributions based on 
QCD sum rules and summarize. 

Our work encloses a number of appendices. Appendix A gives the Fierz transformation 
rules used to derive the isospin relations. Appendix B summarizes the determination of 
nonperturbative parameters in the QCD sum rule approach. Appendix C contains a "hand- 
book" of nucleon distribution amplitudes where we collect all expressions needed in practical 
calculations. 



2 General classification 
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2.1 Lorentz structure 

The notion of hadron distribution amplitudes in general refers to hadron-to- vacuum matrix 
elements of nonlocal operators built of quark and gluon fields at light-like separations. In 
this paper we will deal with the three-quark matrix element 

(0\ e'^u^^z) [a 1 z,a z} i , ti u^(a 2 z) [a 2 z,a z] jlj d^ (a 3 z) [a 3 z, a z} k , tk \P(P,X)) (2.1) 

where \P(P, A)) denotes the proton state with momentum P, P 2 = M 2 and helicity A. u, d 
are the quark-field operators. The Greek letters a, 0, 7 stand for Dirac indices, the Latin 
letters k refer to color, z is an arbitrary light-like vector, z 2 = 0, the at are real numbers. 
The gauge-factors [x, y] are defined as 



[x, y] = P exp 
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ig / &t{x-y)^{tx + {l-t)y) 



(2.2) 



and render the matrix element in ( 2.1|) gauge- invariant. To simplify the notation we will not 
write the gauge-factors explicitely in what follows but imply that they are always present. 

Taking into account Lorentz covariance, spin and parity of the nucleon, the most general 
decomposition of the matrix element in Eq. ( |2.1| ) involves 24 invariant functions: 



A(0\e ijk ui(a 1 z)uj 3 (a 2 z)d k y (a 3 z) \P) = 

= SiMCap ( 75 A0 7 + S 2 M 2 C aP (/ 75 A0 7 + V X M ( 7s C) Q/3 iV 7 + V 2 M 2 ( l5 C) a/3 (^JV) 7 
+Vi (PC) ap ( 75 iV) 7 + V 2 M (fC) ap (/ 75 iV) 7 + V 3 M { ltx C) ap (r 75 iV) 7 
+V4M 2 (JC) ap ( 75 iV) 7 + V 5 M 2 (^C) a p (ia^z^N^ + V 6 M 3 yC) afj (/ 75 iV) 7 
+A! (P 75 C) Q/3 iV 7 + A 2 M (,Pj 5 C) a0 {jN) n + A 3 M { W C) ap (yW) 7 
+A 4 M 2 (/ l5 C) ap N, + A 5 M 2 (j^ 5 C) a/3 {ia^z v N) 7 + A 6 M 3 (fr 5 C) a/3 {/N\ 
+Tt {P»ia, v C) ap (r 75 A0 7 + T 2 M {z»P v io» v C) a p ( 75 iV) 7 + T 3 M (a^C)^ {a^ B N\ 
+%M (P u <T^C) a p {a^z el5 N\ + T 5 M 2 {z v ia, v C) ap ( 7 "75#) 7 
+T 6 M 2 {z^P v ia^C) ap (/ 75 iV) 7 + T 7 M 2 (a^C)^ / 7s iV) 7 

+T 8 M 3 (z"<vC)^ K% 7s iV) 7 , (2.3) 

where iV 7 is the nucleon spinor, C the charge conjugation matrix and a^ u = | [ 7/U , 7l/ ] . The fac- 
tor 4 on the l.h.s. is introduced for later convenience. Each of the 24 functions Si,Vi,Ai, V;, % 
depends on the scalar product P ■ z. 

The invariant functions in Eq. ( |2.3| ) do not have a definite twist yet. For the twist 
classification, it is convenient to go over to the infinite momentum frame. To this end we 
introduce the second light-like vector 

1 M 2 „ . , 
p lt = P lt --z ll , p 2 = 0, 2.4 

2 p ■ z 
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so that P — > p if the nucleoli mass can be neglected M — > 0. Assume for a moment that 
the nucleon moves in the positive e z direction, then p + and z~ are the only nonvanishing 
components of p and z, respectively. The infinite momentum frame can be visualized as the 
limit p + ~ Q — > oo with fixed P ■ z = p • z ~ 1 where Q is the large scale in the process. 
Expanding the matrix element in powers of l/p + introduces the power counting in Q. In 
this language twist counts the suppression in powers of p + . Similarly, the nucleon spinor 
N 7 (P, A) has to be decomposed in "large" and "small" components as 

N~ ( (P, A) = ^— W+ M N,(P, A) = N+(P, A) + N-(P, A) , (2.5) 

where we have introduced two projection operators 

A+ = ^- , A" = -£- (2.6) 
2p-z 2p-z ( J 

that project onto the "plus" and "minus" components of the spinor. Note the useful relations 

p'N(P) = MN + (P), jN(P) = j^N-(P) (2.7) 

that follow readily from the Dirac equation P N(P) = MN(P). 

Using the explicit expressions for N(P) it is easy to see that A + iV = ~ ^fp& while 
A~A r = N~ ~ l/y/p+. To give an example of how such power counting works, decompose 
the Lorentz structure in front of V\ in Eq. (|2.3|) in terms of light-cone vectors 



M 2 M 2 
+ 2j- z i^ N+ ), + (^ iV ") 7 • ( 2 - 8 ) 

The first structure on the r.h.s. is of order (p + ) 3//2 , the second of order (j> + ) 1 ^ 2 , the third of 
order [j> + )~ 1 / 2 and the fourth {p + )~ 3 ^ 2 , respectively. These contributions are interpreted as 
twist-3, twist-4, twist-5 and twist-6, respectively, and it follows that the invariant function 
Vi contributes to all twists starting from the leading one. Such an effect is familiar from 
deep inelastic scattering, where the twist-3 structure function g2(x, Q 2 ) receives the so-called 
Wandzura-Wilczek contribution related to the leading-twist structure function g\(x, Q 2 ). 

The twist classification based on counting of powers of l/p + is mathematically simi- 
lar to the light-cone quantisation approach of ||15|| . In this language, one decomposes the 



quark fields contained in the matrix element Eq. (|2.3j) in 'plus' and 'minus' components 



q = q + + q~ in the same manner as done above with the nucleon spinor Eq. (|2.5| ). The 
leading twist amplitude is identified as the one containing three 'plus' quark fields while 
each 'minus' component introduces one additional unit of twist. Up to possible complica- 
tions due to isospin, one expects, therefore, to find eight independent three quark nucleon 
distribution amplitudes: One corresponding to the twist-3 operator {u + u + d + ), three related 
to the possible twist-4 operators (u + u + d~), (u + u~d + ), (u~u + d + ), three more amplitudes of 
twist-5 of the type (u~u~d + ), (u~u + d~), (u + u~d~) and one amplitude of twist-6 having the 
structure (u~u~d~). 
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twist-3 


twist-4 


twist-5 


twist-6 


vector 


Vi 


V 2 , v 3 


Vk, v 5 


v 6 


pseudo- vector 


Ax 


A 2 , A 3 




A 6 


tensor 


T x 


T 2 , T 3 ,T 7 


T 4 , T 5 , T 8 




scalar 




Si 


S 2 




pseudo-scalar 




Pi 


P2 





Table 1: Twist classification of the distribution amplitudes in Eq. ( |2~§| ) 



Either way, distribution amplitudes of definite twist correspond to the decomposition 
of Eq. ( p.3|) in different light-cone components. After a simple algebra, we arrive at the 
following definition of light-cone nucleon distribution amplitudes: 

4 (0| e^u^axz^a^d^z) \P) = 

= S 1 MC aP ( 75 iV + ) 7 + S 2 MC aP ( 75 iV-) 7 + PiM ( 75 C) a/3 N+ + P 2 M ( l5 C) ap N~ 
+Vi WC)*p (7 5 iV + ) 7 + V 2 WC) ap ( 75 iV-) 7 + yM ( 7± C) Q/3 ( 7 ± 75 iV + ) 7 

v \ r~ m ' 

+Y M (7 x 7 5 iv-) 7 + v 5 — yc) aP ( 75 iv + ) 7 + — v, yc) aP ( 75 iv-) 7 

+Ax (jt l5 C) ap N+ + A 2 ^C) ap N- + (7 ±l5 C) af3 ( 7 ^iV+) 7 

+ Y M (^75^) Q , (7 ± iV~) 7 + A 5 — y l5 C) ap N+ + —As ^C) ap N~ 

+Tt (*cr ±p C) ap ( 7 ± 75 iV + ) 7 + T 2 (i* ±p C) a p ( 7 ± 75 iV-) 7 + T 3 ^ (i^C)^ ( 75 iV+) 

M M 2 M 2 

+T 4 - (.V.pC)^ ( 75 iV-) 7 + T 5 — (i^C) a/3 ( 7 ± 75 iV + ) 7 + — T 6 (i<7 ±z C) Q/3 (tH^T), 



(^ ±±, 7 5 iV + ) + {a L1 _,C) ap (a^'^N 



(2.9) 



where an obvious notation a pz = a^p^Zv etc. is used as a shorthand and _L stands for the 
projection transverse to z,p, e.g. 7 _ l7 - l = ^g^Y with 5^ = ^ - (p^ + z^p v )/pz. 

By power counting we identify three twist-3 distribution amplitudes Vx,Ax,Tx, nine twist- 
4 and twist-5, respectively, and three twist-6 distributions, see Table |I[ Each distribution 
amplitude F = Vi, Ai,Ti, S iy Pi can be represented as 



X ; 



(2.10) 
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where the functions F(xi) depend on the dimensionless variables Xi, < Xi < 1, 'Yli x i = 1 
which correspond to the longitudinal momentum fractions carried by the quarks inside the 
nucleon. The integration measure is defined as 



Jl)x = J dxidx 2 dx 3 5{x\ + x 2 + x 3 — 1) . (2.11) 



Comparision of the expansions in Eq. (|2.3| ) and Eq. ( |2.9| ) then leads to the expressions for 
the invariant functions Vi, Ai, %, Si, Vi in terms of the distribution amplitudes. For scalar 
and pseudo-scalar distributions we have: 



Si = Si, 2p ■ z S 2 = Si - S 2 , 

Vi = Pi, 2p-zV 2 = P 2 -Pi, 



for vector distributions: 



(2.12) 



Vi = Vi, 2p- zV 2 =V 1 -V 2 -V 3 , 

2V 3 = V 3 , Ap ■ zV 4 = -2Vi + V 3 + V 4 + 2V 5 , (2.13) 

Ap ■ z V 5 = V 4 - V 3 , ( 2p ■ z ) 2 V 6 = - Vi + V 2 + V 3 + V 4 + V 5 - V 6 , 
for axial vector distributions: 

A X = A U 2p- zA 2 = -Ai + A 2 - A 3 , 

2A 3 = A 3 , Ap- zA A = -2Ai-A 3 -A 4 + 2A 5 , (2. 14) 

Ap ■ zA 5 =A 3 -A 4 , (2p ■ z) 2 A & = Ai - A 2 + A 3 + A 4 - A 5 + A 6 , 

and, finally, for tensor distributions: 

Ti = Ti, 2p- zT 2 = Ti+T 2 - 2T 3 , 

2T 3 = T 7 , 2p- z% = Ti-T 2 - 2T 7 , 

2p ■ z% = -Ti + T 5 + 2T 8 , (2p ■ zf% = 2T 2 - 2T 3 - 2T 4 + 2T 5 + 2T 7 + 2T 8 
Ap ■ zT 7 = T 7 - T 8 , [2p ■ z) 2 % = -Ti + T 2 + T 5 - T 6 + 2T 7 + 2T 8 . 



(2.15) 



2.2 Symmetry properties 

Not all of the 24 distribution amplitudes in Eq. (|2.9| ) are independent. First of all, each 
distribution amplitude itself has definite symmetry properties. The identity of the two u- 
quarks in the proton together with symmetry properties of quark operators and 7-matrices 
implies that the vector and tensor distribution amplitudes are symmetric, whereas scalar, 
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pseudoscalar and axial-vector distributions are antisymmetric under the interchange of the 
first two arguments: 

^(1,2,3) = ^(2,1,3), 7X1,2,3)= ^(2,1,3), (2.16) 

$(1, 2, 3) = -5,(2, 1, 3) , P(l, 2, 3) = -P 4 (2, 1, 3) , ^(1, 2, 3) = -A t (2, 1, 3) . 

Similar relations hold for the "calligraphic" structures in Eq. (|2.3| ). 

In addition, the full matrix element in Eq. ( |2.9| ) has to fulfill the symmetry relation 

(0|^ fc <(l)^(2)^(3) |P> + (0|^ fc <(l)K(3)4( 2 ) \ p ) + (°i ^X^W^ 1 ) \ p ) = °- 

(2.17) 

that follows from the condition that the nucleon state has isospin 1/2: 

(t 2 - ~(~ + 1)) (0|e^(l)«J(2)d*(3) \P) = 0, (2.18) 

where 

T 2 = - (T+T_ + T_T+) + T 3 2 (2.19) 

and T± are the usual isospin step-up and step-down operators. Applying the set of Fierz 
transformations detailed in Appendix A, we end up with the following six equations: 





27\i 


;i,2,3) = 


[Vi 


-^1(1,3, 2) + [V 1 -A 1 }(2, 3,1), 


[T. + Tj + S, 


~ p i] 1 


;i,2,3) = 


[v 3 


-A 3 ](3,1,2) + [F 2 -A 2 ](2,3,1), 




2T 2 i 


;i,2,3) = 


[T 3 


- T 7 + S 1 + Px](3, 1, 2) + [T 3 - T 7 + 5i + PJ (3, 2, 1) , 


[T 4 + T 8 + S 2 


-P 2 ]i 


;i,2,3) = 


[V 4 


-A 4 ](3,l,2) + [y 5 -A 5 ](2,3,l), 




2T 5 i 


;i,2,3) = 


[T* 


- T 8 + S 2 + P 2 ](3, 1, 2) + [T 4 - T 8 + S 2 + P 2 ](3, 2, 1) , 




2T 6 i 


;i,2,3) = 


[V 6 


-A 6 ](l,3,2) + [y 6 -A 6 ](2,3,l). (2.20) 



The first relation in Eq. (|2.20|) is familiar from (TT| and expresses the tensor nucleon distri- 
bution amplitude of the leading twist in terms of the vector and axial vector distributions. 
Since the latter have different symmetry, they can be combined together to define the single 
independent leading twist-3 nucleon distribution amplitude 



<$> 3 (xi,x 2 ,x 3 ) = [V! - A 1 ](x 1 ,x 2 ,x 3 ) (2.21) 

which is well known and received a lot of attention in the past. The rest of the relations 
in Eq. (|2.20| ) are new and constrain the number of independent distribution amplitudes of 
higher twist. 

In particular, nine distributions amplitudes of twist-4 (cf. Table 1) can be reduced to 
three independent distributions that we choose to be 

$4(xi,x 2 ,x 3 ) = [V 2 -A 2 ] (x 1 ,x 2 , x 3 ) , 

^4(xi,x 2 ,x 3 ) = [V 3 - A 3 ](xi,x 2 ,x 3 ) , 

E 4 ( Xl ,x 2 ,x 3 ) = [T 3 -T 7 + S 1 + P 1 ](x 1 ,x 2 ,x 3 ). (2.22) 
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Lorentz-structure 


Light-cone projection 


nomenclature 


twist-3 


{Cf)®i 


Uj * u d -j- 


$ 3 (^) = [Vi - Ai] { Xi ) 


twist-4 






M x d = l y 2 - M (xt) 




(cy 7± ^) ®tV 




*4(a:<) = [V3 - A 3 ] (xA 






UL -j- %L -|- d 


H 4 (x i ) = [T 3 -T 7 + S 1 + P x \ ( Xi ) 


twist-5 




j i 


$5(2*) = [^5-^5] (Xi) 








*s(a:<) = [Vi -A 4 ] (xi) 








E 5 (x t ) = [T 4 -T 8 + S 2 + P 2 ] ( Xi ) 


twist-6 






M x i) = W* - M 



Table 2: Eight independent nucleon distribution amplitudes that enter the expansion in Eq. ( |2.9|) . 



It is easy to check that all nine distributions V 2 , A 2 , T 2 , V 3 , A3, T 3 , T 7 , Si, P\ can be restored 
from the three above, taking into account the isospin relations in the second and the third 
lines in Eq. flTSCP . 

Similarly, we introduce three independent twist-5 distribution amplitudes 



$ 5 (xi,x 2 ,x 3 ) = [V5 - A 5 ](x 1 ,x 2 ,x 3 ) , 
#5(^1, £2,2:3) = [V4 - A 4 ](x 1 ,x 2 ,x 3 ) , 
55(^1,^2,^3) = [T 4 -T 8 + S 2 + P 2 ](x 1 ,x 2 ,x 3 ) . 

Finally, one distribution amplitude exists of twist-6: 

$ 6 (xi,x 2 ,x 3 ) = [V e - A 6 ](xi,x 2 ,x 3 ) . 



[2.23) 



[2.24) 



2.3 Representation in terms of chiral fields 

The physics interpretation of the distribution amplitudes $(xj) is most transparent in terms 
of quark fields of definite chirality 



q 



TU) 



1 



i±7 5 )g. 



(2.25) 



Projection on the state where the spins of the two up-quarks are antiparallel v/u^ singles out 
vector and axial vector amplitudes, while parallel spins v)u^ or correspond to scalar, 
pseudoscalar and tensor structures. We are lead to the classification of the distribution 
amplitudes in terms of the quarks light-cone components and spin projections, summarized 
in Table £| The leading twist-3 distribution amplitude can be defined as (cf. ||16||): 

(Qle^fulMCMjM^) id\{a 3 z) \P) = ~pz ^ jvx e~ ipz ^ **« $ 3 ( Xi ) . (2.26) 



S 



Twist-4 distributions allow for the following representation: 

(0 1 e ljk (uj (a lZ )C tu\ (a 2 z) ) i>d\ (a 3 z) \P) = - ^pz ^ J Vx e' tpz ^ x * a < $ 4 (x t ) , 

(0| e ijk (uj{a lZ )C /yjj^oa-z)) j x ^d l k {a 3 z) \P) = -pzM/N^ jVxe~ ipz ^ Xiai * 4 (»i) , 

(0\e» k ful(a lZ )C ^u](a 2 z)j k (a 3 z)\P) = ^/c.\//.V' Jvxe~ ipz ^ x ^ S 4 (x i )(2.27) 
and, similar, for twist-5: 

(0| £^(4(0x^(7^(02^) /4(a 3 ^) |P) = -^MViV T jvxe~ ipz ^ x ^ $ 5 (^) , 

(0|£^(wT(a 1 z)C^7 ± /wj(a 2 z))7V4( a 3^)| J P> = -pzMp'N^ J Vxe~ ipz ^ Xiai # 5 (^) , 

(0|e y *(uJ( ai z)C7^ut(a 2 z))j«J(a3z)|P) = -pzMtfrf Jvxe~ ipz ^ x ^ E 5 (xi) (2.28) 
Finally, the twist-6 distribution amplitude is written as 

r iifc (uI(aiz)Ctyij(a 2 s)) j*4(a 3 z) |P> = -^M 2 ^V T Jvxe~ %pz ^ x ^ $ 6 (x 4 ) . (2.29) 



In the rest of the paper we study this set of distribution amplitudes in some more de- 
tail. Final expressions for all the 24 distribution amplitudes in Eq. ( |2.9| ) are collected in 
Appendix C. 

3 Conformal expansion 

The conformal expansion of light-cone distribution amplitudes is the field-theoretic analogue 
to the partial wave expansion in quantum mechanics. The idea, in both cases, is to use the 
symmetry of the problem to introduce a set of separated coordinates. In quantum mechanics, 
spherical symmetry of the potential allows to separate dependence on radial coordinates 
from angular ones. All angular dependence is included in spherical harmonics which form 
an irreducible representation of the symmetry group 0(3), while the dependence on radial 
coordinates is governed by a one-dimensional Schrodinger equation. 



In the same spirit conformal symmetry |L7j of the QCD Lagrangian can be used to 



study the distribution amplitudes as it allows to separate longitudinal degrees of freedom 



from transverse ones |L8L [L9L [21], H |22], 16 . The dependence on longitudinal momentum 



fractions is taken into account by a set of orthogonal polynomials that form an irreducible 
representation of the collinear subgroup SL(2, R) of the conformal group which describes 
Mobius transformations on the light-cone. Transverse coordinates are replaced by the renor- 
malization scale, the dependence on which is governed by the renormalization group. Since 
the renormalization group equations to leading logarithmic accuracy are driven by tree-level 
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counterterms, they have the conformal symmetry. As a consequence, components in the 
distribution amplitudes with different conformal spin, dubbed conformal partial waves, do 
not mix under renormalization to this accuracy. 
Conformal spin of the quark is defined as 

3 = \{l + s) (3.1) 

where I = 3/2 is the canonical dimension of the quark field and s = ±1/2 is the quark 
spin projection on the light-cone. The spin projection operators are in fact the same as 
used to separate the "plus" and "minus" components of a spinor in Eq. (|2.6|) . The "plus" 
component of the quark field q + corresponds to s = 1/2 and, therefore, j = 1, while the 
"minus" component q~ has s = —1/2 and j = 1/2. For multiquark states, we are left with 
the classical problem of summation of spins, with the difference that the group is in our 
case non-compact. The distribution amplitude corresponding to the lowest conformal spin 
jmin = ji + 32 + 33 of the three-quark system is equal to 



Contributions with higher conformal spin j = j min + n,n = 1,2, .. . are given by $ as mul- 
tiplied by polynomials which are orthogonal over the weight function Eq. (|3.2| ). A suitable 
orthonormal basis of such "conformal polynomials" has been constructed in |16| . 



In this section we consider the conformal expansion of nucleon distribution amplitudes 
taking into account contributions of leading and next-to-leading conformal spin (i.e. "S" 
and "P"- waves). The expansions can easily be extended to arbitrary spin. One reason why 
we do not present complete expansions in this paper is that we do not have the tools to 
estimate the corresponding additional parameters. The second reason is that for yet higher 
spins one has to take into account contributions of four-particle distributions involving an 
extra gluon that we do not consider here. 

3.1 Leading twist-3 distribution amplitude 

At leading twist, there is only one independent distribution amplitude corresponding to the 
light-cone projection (C f) ® / and therefore involving three "plus" quark fields , see Table 2. 
The conformal expansion then reads 

$ 3 (xi,fj) = 120xix 2 x 3 [<fM + 0j(/i)(xi - x 2 ) + 03 (/i)(l - 3x 3 ) + . . .] . (3.3) 

Here /i is the renormalization scale. For the discussion of the shape of the distribution it is 
often convenient to rewrite Eq. ( |3.3j ) to factor out the overall normalization: 



<$> 3 {xi,fj) = 12Oxix 2 x 3 03(/i) l + 3 ( / u)(xi-X2) + 03 f (/i)(l-3x 3 ) + ... . (3.4) 



The relation between (ftf and 4>f is obvious. 
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The statement of conformal symmetry is that the coefficients ± (/i) do not mix with 
under renormalization since they have different conformal spin: j = 4 and j = 3, 



respectively. By an explicit calculation one obtains [23, U% ER Eq 



0°(^ 2 ) = L*'*><jP 3 M , L ee ^ , (3.5) 
where 6=11 — 2/3n/, and 

»M = ~ (3^ 20/% + £ 8/3b ) + J (^ 20/9fe - ^ 8/3b ) & W , 

03 W = ^ (L 2 ^ - L 8 / 3 *) g^) + 1 (L 20 ^ + 3L 8 / 36 ) ^) . (3.6) 

Numerical estimates for the coefficients are available from QCD sum rules [|IT| :Q 

0°Gu = 1 GeV) = (5.3 ± 0.5) • 10~ 3 GeV 2 , 
0+(/i = 1 GeV) = 1.1 ±0.3, 

0- (^ = 1 GeV) = 4.0 ±1.5. (3.7) 



More sophisticated models suggested in [11] involve in addition contributions of second 



order polynomials related to the operators with conformal spin-5. The estimates of the 
corresponding coefficients are, however, less reliable and have large errors. 

3.2 Higher-twist distribution amplitudes 

The conformal expansion of the higher- twist distribution amplitudes defined in Sect. 2 is 
equally straightforward. With the help of the general expression in Eq. (|3.2|) one obtains for 

$4(^1) = 24XiX 2 [0° + 07 (Xi - X 2 ) + 04 (1 - 5x 3 )] , 

ty 4 (xi) = 24xix 3 [V>2 + ^{x! - x 3 ) + ^ 4 + (l - 5x 2 )] , 
S 4 (x i ) = 24x 2 x 3 + Zl(x 2 - x 3 ) + £+(1 - 5x x )] > (3-8) 



twist-4: 



for twist-5: 



and for twist- 6: 



$5 fai) = 6X3 [0° + 05 (Xi - X 2 ) + 0^(1 - 2X 3 

V 5 ( Xi ) = 6x 2 [iP° 5 + ^ 5 -(x! - x 3 ) + ^ 5 + (l - 2x 2 ) 
Zs(xi) = 6x1 fe° + e 5 "(^2 - x 3 ) + - 2xi)] , (3.9) 



%(x t ) = 2 [0° + 06 (xx - x 2 ) + 0+(l - 3x 3 )] • (3.10) 

'In notations of = /jy. The given numbers correspond to the last reference in [ |TT| . 
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At this point the expansion introduces 21 new parameters. Our next task is to find out how 
many of the parameters actually are independent and which are connected by equations of 
motion. 

The normalisation of all distribution amplitudes and, therefore, the asymptotic wave 
functions are determined by matrix elements of a local three-quark operator without deriva- 
tives. The Lorenz decompositon of a local three quark matrix element is much simpler 
compared to the general parametrisation in Eq. (|2.3|) and involves only four structures: 



4 (0| e «X(°)4(°)d*(°) \N, P) = V? (PC) af3 ( 75 A0 7 + V 3 °M ( 7 ,C) Q/3 ( 7 " 75 A0 7 

T° {P v io» v C) a(i ( 7 "7 5 A0 7 + %°M (<vC)o0 K*7 5 A0 7 • (3-11) 

From isospin constraints it follows that in addition V° = . Thus there exist only three 
independent constants. 

Remarkably, these three parameters are well known and can be obtained from existing 
estimates of the following three matrix elements: 

(0| e^ k [u\0)C MO)] 75 #*(0) \P) = hpz jtN(P) , 
(0\e^ k [^(O)CV^'(O)] 7s7^ fc (0) \P) = AiMJV(P) , 
(0| e ijk [^(0)0(7^(0)] l5 a^d k (0) \P) = X 2 MN(P) . (3.12) 

The parameter /n = V° enters already at the level of twist-3 and determines the normal- 
isation of the leading twist distribution amplitude Eq. (|3.7| ). The other two parameters 
Ai = (V° — 4V° ) and A 2 = 6(V° — 4Zj°) correspond to the nucleon coupling to the two possi- 
ble independent nucleon interpolating fields that are widely used in calculations of dynamical 
characteristics of the nucleon in the QCD sum rule approach. The operator corresponding 



to Ai was introduced in |T2] while the one corresponding to A2 was advertised in |TJ 



The QCD sum rules summarized in App. yield the following estimates: 

AiGu = 1 GeV) = -(2.7 ± 0.9) ■ 10~ 2 GeV 2 , 

A 2 (/i = 1 GeV) = (5.1 ± 1.9) • KT 2 GeV 2 , (3.13) 
Anomalous dimensions are the same for both currents p7| 



X 1 (ji 2 ) = L 2/ %(/i!) , A 2 (^ 2 ) = L 2 /%Gu 2 ) . (3.14) 

Note that one overall sign in the determination of vacuum-to-nucleon matrix elements is 
arbitrary as it corresponds to arbitrary (unphysical) overall phase of the nucleon wave func- 
tion. The relative signs between the couplings /at, Ai, A 2 , etc. are, however, well defined and 
can be determined from suitable nondiagonal correlation functions, see App. |TJ for details. 
We choose /n to be real and positive; then, it turns out that Ai is negative and A 2 positive. 



The fact that Ai and A 2 have opposite signs is known from |14l , the negative relative sign 
between Ai and /n is a new result. 
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The coefficients in Eqs. (|3] 
conformal spin then read 



15[), ( |3.10D corresponding to the operators of leading 



CO _ cO 
S4 _ S5 



1 



1 

2 



6 



-A, 



$ = ^°=2 (/" - Ai) 



Note that the normalization of twist-3 and twist-6 distribution amplitudes are equal, and 
similarly for twist-4 and twist-5 distributions. The constants and 1P45 involve both 
the leading twist matrix element Jn Eq. ( |3.12| ) and the higher- twist contribution ~ Ai. 
The former is analogous to the Wandzura-Wilczek-type contribution to the higher-twist 
distribution amplitudes studied in || for the p-meson, and the latter defines the "genuine" 
higher- twist correction. Note that Ai, A2 /at that is in agreement with our expectation 
(see Introduction) that the matrix elements of higher-twist three-quark operators are large. 

The remaining contributions of next-to-leading conformal spin are related to operators 
with one derivative. In much the same way as before and as elaborated in App. || the 
14 unknown parameters are reduced by isospin symmetry and equations of motion to five 
dimensionless parameters Vf, A™, ff, / 2 d , /" which we define in the following way: 

(0| (u(0)C7u(0)) 75 j(izl5d)(0) \P) 



u(0)Cfy 5 iz Du(0)j /d(0) 

(0| («(0)C 7M tt(0)) 7 57 M j(izl}d)(0) 
(u°(0)C<v«(0)) 75(7^ ^tz^d)(0) 
'u(0)C W iz D u{0)) ^ ^d(O) 



P) = 


hVfipzf jN(P) , 


P) = 


-hAlipzf tN{F) , 


P) = 


\xft(pz)M JN(P) 


P) = 


-\ 2 ffoz)M fN{P) 


P) = 


-\ 1 f?(pz)M Z'N(P) 



(3.15) 



Here, we have used the shorthand notation for the left-right derivative iz- D= iz ■ (D — D) 
and for brevity omitted colour indices. The first two matrix elements and A\ are leading 
twist-3 and were estimated in flTT]p} 



1 GeV) 
1 GeV) 



0.23 ±0.03, 
0.38 ±0.15. 



(3.16) 



We have used these estimates above in Eq. ( |3.7| ): 4>^ = ^A\ anc ^ 03 = 2O- ~ ^Vf)- The 
other three parameters are genuine higher-twist and are estimated in App. |Bj using QCD 
sum rules. We obtain: 



1 GeV) 
1 GeV) 
1 GeV) 



0.6 ±0.2, 
0.15 ±0.06, 
0.22 ±0.15. 



(3.17) 



Uln notation of (TJ Vf = V^^.A^ = 2A^ ' 1 



.0) 
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The remaining coefficients expressed in terms of the above parameters read for twist-4: 
C = \ (Ai(l - 2ff - Aft) + f N {2Ai - 1)) , 

= ~ (Ai(3 - lOtf) - f N {lW? - 3)) , 

ft = ~\ (Ax (2 - 7tf + ) + MA£ + 31? - 2)) , 

ft = ~\ (Ai(-2 + hfl + 5 ft) + f N (2 + 5A« - SVf)) , 

C = ^A 2 (4 - 15/ 2 d ) , 

£ 4 + = ^A 2 (4 - 15/ 2 d ) . (3.18) 



for twist-5: 



|(A 1 (/ 1 d -/r) + /iv(2^-l)) , 
-^(X l (Aft-l) + f N (3 + AV 1 d )) , 



and for twist-6: 



ft = \ (Ai(/i d - fi) + /at(2 - A? - 3V*)) , 
ft = — (Ai(-1 + + ft) + f N (5 + 2A\ - 2V?)) , 

£+ = ^A 2 (2-9/ 2 rf ), (3.19) 
06= ^(A 1 (l-4/ 1 d -2/ 1 «) + /^(l + 4Ar)) , 

6 + = (Ai(l - 2ff) + /^V? - 1)) . (3.20) 

With these relations, the construction of higher-twist nucleon distribution amplitudes is 
complete to our accuracy. Note that numerical values of the parameters given in Eqs. (|3.16|) , 
( [3.13] ) are strongly correlated within the QCD sum rule approach, and the given errors should 
not be added together in the sums fl3.18|) -( |3.20D . We expect that the accuracy of the QCD 
sum rule calculation of all coefficients in the distribution amplitudes is of the order of 50%. 
In table Tab. (0) we have collected the numerical values for the expansion parameters. 

In Fig. [I] we have presented the leading twist distribution amplitude $ 3 (a;j), in Fig. |2] the 
three distribution amplitudes of twist-4 are plotted. Note that the scale on the z-axis is not 
the same for each plot. 
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Figure 1: Twist-3 distribution amplitude $3(3^) 



4 Summary and conclusions 

We have carried out a systematic study of the higher-twist light-cone distribution ampli- 
tudes of the nucleon in QCD and found that a generic three-quark matrix element on the 
light-cone can be parametrized in terms of eight independent nucleon distribution ampli- 
tudes. In particular, we have identified one distribution of leading twist-3, three of twist-4, 
three of twist-5 and one of twist-6. The shape of the distribution amplitudes at asymtoti- 
cally large momentum transfers is found and is dictated by the conformal symmetry of the 
QCD Lagrangian. In order to quantify the corrections, we attempt to expand the distribu- 
tion amplitudes in contributions of local conformal operators and keep contributions of the 
next-to- leading conformal spin ("P" -waves). The coefficients in this expansion present the 
necessary nonperturbative input. Some of them are related by QCD equations of motion, so 
that to our accuracy we end up with eight independent parameters 

From this amount, three parameters determine the leading twist amplitude and are known 



already from the analysis in |nj ; two more — the normalisations of the higher- twist asymp- 
totic distributions — can be directly inferred from QCD sum rule phenomenology of the 
nucleon. The remaining three parameters determine the deviation of the higher-twist distri- 
bution amplitudes from their asymptotic form and are estimated in this work. 

To avoid confusion, we remind that our work does not present the complete analysis 
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Figure 2: Twist-4 distribution amplitudes $4(2;*) in the first line to the left, ^4(2^) in the first 
line to the right, and H(xj) in the second line 

of all existing higher-twist distributions, but of its subset related to three-quark operators 
without extra gluon fields. As explained in the introduction, we believe that the three- 
quark contributions considered in this paper dominate higher-twist corrections because of 
large matrix elements. Indeed, we found that the matrix elements of higher-twist three-quark 
operators are of the same order or larger than those of the leading twist, see Sec. 3. We stress 
that the situation in the nucleon case is quite different to that of higher twist-corrections in 
the meson case studied earlier || |]. The reason is that in the meson case "bad" components 
of quark fields always can be eliminated in favour of "good" components and an additional 
gluon. The resulting matrix elements turn out to be numerical small. In the nucleon case 
quark fields with "minus" projection give rise to genuine three-quark higher twist effects 
that cannot be attributed to higher Fock components in the wave function, involving an 
additional gluon. 

A detailed summary of the distribution amplitudes is presented in Appendix C. In this 
paper we do not consider phenomenological applications, but expect that our results are 
relevant for the studies of a wide range of interesting physical processes. 

One obvious application can be found in the calculation of the nucleon form factors. While 
the magnetic form factor of the proton Gm(Q 2 ) allows for the leading twist description |3j. 
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it is known that the Pauli formfactor ^(Q 2 ) is higher-twist, suppressed by an additional 
power of Q 2 . The description of F 2 therefore calls for the input of higher twist distribution 
amplitudes presented in this paper. 

Twist-4 distribution amplitudes of the nucleon are, probably, the most interesting and in 
general are related to various spin asymmetries in exclusive processes. Indeed, while the total 
helicity is conserved in processes involving leading twist distribution amplitdues |29| , this is 



not the case for non-leading twist contributions. Thus spin-sensitive transition formfactors 
will provide a testing ground for higher twist effects. 

Last but not least, the distribution amplitudes presented in this paper can immediately be 
applied to calculations of exclusive reactions at moderate momentum transfers in the frame- 



work of QCD light-cone sum rules J30|. This approach allows to study soft non-factorisable 



contributions to hard exclusive reactions in a largely model-independent way, see e.g. pi 
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Appendices 

A Fierz transformations 

In this section we give the Fierz transformation rules necessary to derive the isospin con- 
straints Eq. ( |2.20| ) from the symmetry requirement Eq. ( p . 1 7| ) . 

We write the definition of distribution amplitudes in Eq. ( |2.9|) in a shorthand notation 

as 

<o| ^X(^4( a ^K(M \p(p, A)) 

= Si( + s 2 ( s 2)q/3,7 + -Pl(Pl)a/3,7 + ^2(^2) a/3, 7 + Vl ( v l ) a/3,7 ^2(^2)0/3,7 7> ^3(^3)0/3,7 



■iv 4 (f 4 ) a/ 3 i7 + V 5 (v 5 ) a p a + V 6 (v 6 ) al3 „ + A 1 (a 1 ) aPa + A 2 (a 2 ) Q /3, 7 + 2^ 



+—^.4(04)0,3,7 + ^5(05)0/3,7 + ^6(06)«/3,7 + Ti(ti) a j3 a + T-zit-ijap^ + T^{t^) a p^ 

+T 4 (t 4 ) a p n + T 5 (t 5 ) a p n + -T 7 (t 7 ) a p n + -T 8 (t$) a p tl . (A.l, 

The small letters, e.g. {s\) a p^ = {C) a pN 1 stand for the Lorentz structures (TC) a p (T'N) 
etc. By means of the following Fierz transformation 

(m a p (r'iv) 7 = \ \c 7 p (rr'N) a + ( 75 c) 7/3 (r 75 r'iv) a 
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+^r z <&) 7fi (r F'N) a + ± yc) lP (r pT'N) a + ( 7J c) 7/J (r^r N ) a 

-^- z ^c) lP (r ^v'N) a - J- (^ 75 c<) 7/3 (r ^v'N) a - ( 1±l5 c)^ {r 1 ± l5 r'N) a 



^ (a pz C)^ (Ta zp T'N) a + \ {a^C\ p (iV^'r'A^ 



(P 

+^7 (^) 7 /3 fl*N) a + -L. {ia^ z C) lP (iV pT'N)] (A.2) 

all Lorentz structures are brought into the same form, i.e. we apply for the expansion of the 
third contribution in Eq. (|2.17 ) the following transformations of twist-3 structures 

(oi) 7 Aa = ^ + ai - ti) aA7 

(ti) 7 /3,a = - (vi + ai) aPl (A.3) 
and similar transformations for twist-4 



a/3,7 



a/3,7 



(si) 7 /3, a = ^ ( s x + pi - 2w 2 - 2a 2 + v 3 -a 3 -t 3 - 2t 2 + ~t 7 

(Pl)"fP,a = ~ f Si + Pi + 2f 2 + 2fl 2 - ^3 + A3 - *3 ~ 2t 2 + ~*7^ 
(^) 7 /3,a = J f-Si + Pi + V 3 + a 3 - t 3 - ^ 7 J 



(02)7/3,0 = ^ ( _s i + V\ ~ v 3 ~ a 3 - t 3 - -t 7 



a/8,7 



0/3,7 



(^3)7/3,0 = 7 (2si - 2p a + 4t> 2 - 4a 2 - 2t 3 - t 7 ) a/3j7 
(03)7/3,0 = 7 (~2si + 2pi + Av 2 - 4a 2 + 2t 3 + t 7 ) a 



v ^yi -r ^2 -r ^"3 t "7; a/ 3 )7 

1 

4 



(t 2 ) 7 /3, a = 7 (-2si - 2 Pl - 2t 3 + t 7 ) ap 



(t 3 )-y(3,a = ~ ("Si -Pi- 2^2 - 2a 2 - V 3 + CL 3 + t 3 - 2t 2 ~ ~t 7 



a/3,7 



(t 7 )yf>, a = \ {2s 1 + 2 Pl - Av 2 - 4a 2 - 2v 3 + 2a 3 - 2t 3 + At 2 + t 7 ) a/3 , 7 • (A.4) 

Relations for twist-5 and twist-6 are identical up to the substitutions: {v±, ai, ti} — > {v 6 , a 6 , £ 6 } 
and {si,pi,v 2 ,v 3 ,a 2l a 3 ,t 2 ,t 3 ,t 7 } — > {s 2j p 2j v 5 ,v 4j a 5 ,a Aj t 5 ,t 4 ,t 8 }. By a simple substitution, 
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one also gets the necessary Fierz transformations from (. . .) a -y,p -*•(•■ •) a/3,7 required in the 
second contribution in Eq. (|2.17|) . Putting everything together we derive three equations for 
the twist-3 amplitudes from the coefficients Vi,ai,t\: 

= [271(1, 2, 3) - V x (l, 3, 2) + A 1 (l, 3, 2) - V(3, 2, 1) - A 1 (3, 2, 1)] {t x ) aPn 
+ [2V! (1, 2, 3) + 2Vi(l, 3, 2) + 2^(1, 3, 2) - 271(1, 3, 2) 

+V1(3, 2, 1) - ^(3, 2, 1) - 271(3, 2, 1)] (v x ) a ^ 
+ [2^(1, 2, 3) + 2^(1, 3, 2) + 2^(1, 3, 2) + 271(1, 3, 2) 

-V(3, 2, 1) + A 1 (3, 2, 1) - 271(3, 2, 1)] (a^ . (A.5) 

Using the symmetry properties Eq. ( |2.16| ) it is easy to see, however, that all three above 
equations are in fact identical to the one given in Eq. (|2.20|) . In the similar way one can 
derive an overcomplete set of nine twist-4 equations that all can be reduced to the two 
equations relating twist-4 amplitudes in Eq. (|2.20|) . The derivation for the remaining twist-5 
and twist-6 equations is similar. 

B Nonperturbative parameters 

In Sect. 3.2 we have demonstrated that the normalisation of the asymptotic distribution 
amplitudes of all twists involves three independent parameters. In the first part of this 
Appendix we consider leading corrections to the asymptotic distribution amplitudes, related 
to contributions of three-quark conformal one extra derivative, and find that they involve 
five new parameters. 

In the second part, we work out QCD sum rule estimates of the higher-twist matrix 
elements. Partially the results can be adapted from existing sum rule calculations containing 



nucleon currents [121 14, E8|, partially they are new. 



B.l Equations of motion 

The coefficients of "P"-wave contributions in the conformal expansion of distribution am- 
plitudes are given by matrix elements of operators containing one derivative. In case the 
derivative acts on the down quark we have the general decomposition: 

4 <0| ^V a (0)4(0) [iD^f (0) \N, P) 

= VfPx (PC) a(3 ( 75 A0 7 + V 2 °M (PC) ap ( 7 a7 5 A0 7 + VlP x M ( 7m C)^ ( 7 " 75 iV) 7 
+V 4 °M 2 { lx C) ap ( 75 iV) 7 + V 5 °M 2 (YC) ap (^ A75 iV) 7 

+T 1 d P x {P v ia^C) ap (r 75 A0 7 + %°M (P u i<TXvC) af} ( 7s iV) 7 + T 3 d P x M {o, v C) a(i (^ 75 iV) 7 
+T 4 °M (P^C) ap (^ A75 iV) 7 + T 5 °M 2 (ia»xC) af3 ( 7 " 75 iV) 7 

+T 7 °Af 2 (<vC) a/3 (a^ 7A75 iV) 7 (B.l) 
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yielding 11 parameters. 

Acting with one derivative on the up-quarks picks up the structures of 7-matrices that 
are odd under transposition: 
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<{0)iDx u 3 p (p)d*(p)\N,P) 



= S?P x MC af3 ( 75 iV) 7 + S»M 2 C af3 ( 7 a7sA0 7 + 
+V?P X M ( 7s C) a/3 iV 7 + P 2 °M 2 ( 75 C) Q/3 ( 7A iV) 7 + 

+A U ,P X (Pl 5 C) a p N, + A%M (P 75 C) Q/3 ( 7 AiV) 7 + AIP X M {l,l,C) a p ( 7 ^A0 7 

+A° 4 M 2 ( lxl5 C) a(3 N, + A^M 2 (YlsCr? («r M JV) 7 (B.2) 

and introduces 9 parameters. Hence, altogether there are 20 unknown numbers. The equa- 
tions of motions yield, however, a number of constraints: 



(0| £^V(0)CVi J '(0) 7 A [iD x d^} k (0) \N, P) = 



(0| e %3 V(0)C 7 V(0) [iD A rf 7 r (0) |AT, P) = P x (0| e« V(0)C 7A u'(0) d*(0) |iV, P> , 
(0| ^V(0)Ccr a/3 ^(0) 7 A [il^]* (0) |JV, P) = , 
(0| e* ik u i (0)Cia al3 u'(0) [iD^d^] k (0) |iV, P) 

= P' 3 (0| e iik u\G)Ci<T a pv?(G) d k (0) \N, P) - (0| £^' fc 

(0| £^V(0)« 75 <7 a/3 ^'(0) [iD^]* (0) |JV, P) 

= P^ (0| e ij 'V(0)C7 5 ia^M J (0) d k {0) \N, P) 



u(0)CtD a u(0) d k J0)\N,P) 



- u e 



w(0)C 7 5« P a w(0)j <(0)|iV,P), 

(B.3) 



which translate to 



Vf = 4V 2 ° + 2V 3 d . 



-3V 5 ° = Vl + V° 



4 • 



y° - V 3 ° = V? - V 3 d + 4V 4 ° - V 2 C 

= T 3 d + T 5 ° , T 2 ° = -T? + 47"/ + 3T 4 ° , 

T* - 2T 3 d - T 4 ° + 3T 5 ° + 6T 7 ° , 
T x - 2T 3 ° - = 7? - 3T 2 ° - 2T 3 d - T 4 ° , 

-27/ + 2T 4 ° = -2T 3 ° - VI , 2T 3 ° - P 2 ° = 27/ - 2T 4 ° - 6T 7 ° , 



7J° - 27/ - & 



(B.4) 

The two relations in the first line in Eq. ( |B.4| ) follow from the first equation in Eq. ( |B.3| ) , the 
relations in the second line from the second equation, the third line from the third equation, 
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fourth and fifth line from the fourth equation and the last line from the last equation. In a 
similar way one observes that the axial structures have to fulfill the constraints 



Sjk 



Kl 



u(0)C 7 « 7fi D e u(0) d k (0)\N,P) = 



T ijk 



u(0)C 7 A« D e -l e i Dx 75«(0) <(0) \N, P) 



(0\e^ k 



u(0)C- \ a XQl a iD a + j a a Xe iD a \ l5 u(0) 



1 y 



d k {0)\N,P) 



P a (0| e ijk u i (0)C^ s u j (0)d k (0) \N, P)-(0\ ^V(0)C 7 V(0) [iD a d 1 ] k (0) \N, P) 

(B.5) 

where from it follows that 

A\ + At + A° 2 + AAl = , 

Al-A Q 2 = V 2 ° + V° - V d , , 2A° 5 = V 2 ° + V 3 ° - 2 V 5 ° - V 3 d . (B. 6) 

where the relation in the first line is derived from the first equation in Eq. ( p.5|) . Combining 
the constraints from the equations of motion with the eight isospin relations obtainable from 
Eq. (|2.20|) we obtain an overcomplete set of 20 equations. Choosing Vf, Af, V| , A3 and T 3 d 
as the independent parameters and solving these equations, we can express the remaining 
parameters as 



Vl = \{V d 1 -2Vt) 



l(-4V 1 + 3V 1 d + 4V 3 -18V 3 d ) 



V° = — (4 V? - 3V? - 4 V 3 ° + 2 V 3 d ) 
16 

A» = At-^-V! + ^, 



^ = _L i-AA\ -8A% + Vf + 4 V 3 ° - 6 Vf) , A °5 = ^ ( 4V i° + 3 V i + 20 V 3° _ 18 V s) 



1 



(Al + 6 At + 8 T* - V? - 2 V$ - 6 V 3 ° + 12 V 3 d ) , T 5 ° = -T 3 



d 

3 ' 



T 4 ° = - (-A u x + 2 At -8T 3 d + V° - 2 V$ - 2 V 3 ° + 4 V d ) , 
T 7 ° = — (ZA\ - 6 At - 16T 3 ° + 48T 3 d + V? + 6 Vf + 6 V 3 ° - 12 V 3 d ) , 
Si = — (3A1 + 10 At - 16T 3 ° + 16T/ + 9 Vf - 2V? - 10 V^ + 20V 3 d ) 
S« = - (3A± + 10 A% - 8T 3 ° + 16T 3 d + V? -2Vf- 10 V 3 ° + 20 V 3 ) , 
VI = j (Al - 2 At - 8T 3 ° + 16T 3 d - V? + 2 V? + 2 V 3 ° - 4V 3 d ) , 
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Vl = — {-A\ + 2AI + 16T 3 ° - 16T 3 d + 5 V? - 2Vf - 2 V 3 ° + 4 Vf) 
^ = 5 Mi + V? - V?) . 



(B.7) 



Finally, the parameters Vf, A±, Vf, -4. 3 and 7^ can be expressed by the reduced matrix 
elements of the operators defined in Eq. (|3.15|) . We have 



Vf = f N Vf , A\ = f N A\ , XJl = - (Vl - 6V 3 d ) , 

Ai/r = ~ (-2AI - 12AI + Vl + 4V 3 ° - 6V 3 d ) , 
A 2 / 2 d = 2 (-2A% - 87? + Vl + 2V 3 ° - 4V| + 27?) . 



(B.8) 



The first two parameters appear already at the level of twist-3 and are known from ||11|| . The 
others will be calculated in the next subsection. 



B.2 QCD sum rule estimates 

The QCD sum rule estimates for Ai, A 2 are derived from the consideration of the two-point 
correlation function 

1 J d 4 xe^^O\r ]l (x)m\0) = \X i \ 2 M 2 ^^l + ..., (B.9) 

where (0|t?i(0) \P) = MXxN(P) and (0\r) 2 (0) \P) = MX 2 N{P)^\ are the local three-quark 
operators defined in Eq. ( [3 . 1 2j ) . The dots refer to excited states and the continuum contri- 
bution. Taking into account vacuum condensates up to dimension-8 one obtains the sum 



rules 14, 28 



2(2vr) 4 M 2 |A 1 | 2 = exp(M 2 /Ml) ^M B E 3 (s /M B ) + -M^so/Ml) + - U - j± 

(BIO) 

and 

2 (2tt) 4 M 2 M = exp(M 2 /M 2 ) {m b E 3 (s /M b ) + ^M^so / >M 2 B )\ , (B.ll) 

where 

n-i 1 / k 

En{s»lMl) = l-e^l M ^h{w) ■ (R12) 

fc=0 ^ B ' 



f Note that our normalisation of Ai differs by a factor M from the standard one 
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Here and below Mb is the Borel parameter. We use the Borel window lGeV 2 < M B < 
2GeV 2 , with the continuum threshold t/so ~ 1.5 GeV and values of the condensates nor- 
malized at /i 2 = 1 GeV 2 @: 

a = -(2n) 2 (qq) ~ 0.55 GeV 3 , 
b= (27r) 2 (— G 2 ) ~ 0.47 GeV 4 , 



7T 



m 2 = -0.8 GeV 2 . (B.13) 

w) 

With these inputs, we find |A X | = 0.027 ± 0.009 GeV 2 , |A 2 | = 0.051 ± 0.019 GeV 2 . Note 
that the above sum rules only fix the absolute value of the parameters Ax and A 2 . Relative 
phases between different nucleon-to-vacuum matrix elements can be computed by investi- 
gating suitable non-diagonal correlation functions. To determine the relative sign between 
Jn and Ai we consider the correlation function 

i J d A xe^ x (0\e^ k [u\x)C^{x)] 7 5^ fc (x) ^(0) |0) = /^M/ ^j^ + . . . .(B.14) 

Taking the ratio of the corresponding sum rule and the sum rule Eq. ( |B.10| ) we obtain 

h = _l (2M 2 B E 1 (s /M 2 )-m 2 )aM 

Al 3 M B E 3 (s /M 2 ) + h -Ml El {s,/Ml ) + £ U - 

which is real and negative. We have chosen f N to be positive according to the standard 
choice in []nj and Ai to be negative. Note also that the above sum rule leads to a value /at = 
5.6 ■ 10~ 3 GeV 2 consistent with standard estimates. A similar calculation of the nondiagonal 
correlation function 

i J d 4 ze^ (0| TftOr^O) + TfeOsOTft (0) |0) = (AxA* + A 2 A*) M 2 ^^ 2 +... (B.16) 



was performed in [14] . The resulting sum rule 



- 2(2vr) 4 M d (AiA; + A 2 A^) = 2ae M B 12M B E 2 (s /M 2 B ) - Qm^E^ / M 2 B ) + n _ f > 

\ " 2l~.\lf, 

(B.17) 

shows that the relative sign between A 2 and Ai is negative. 

To estimate the remaining parameters we consider the nondiagonal correlation functions 
involving the last three higher-twist operators in Eq. (|3.15| ) with either 771(0) or 772(0): 

/— f d \ A I 2 M 2 a 7 d 

d 4 xe^ (0| e^ k (^(x)C7^(x)) lbl ^iz^d) k {x) ^(0) |0) = ^ 2 _ g2 + • • • > 
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i I d 4 xe^ (0|^ fc {u{x)C llxl5 izDu{x)Y ^d\x) ^(0) |0) = f ^ X ^ 9 2 — + ■ ■ ■ 



i I d 4 xe^ fc (0| {u\x)Ca, v u>(x)) l5 a^(izDd)\x) fj 2 {0) |0) = H 1 2 ' \ g + 



M 2 -q 2 

\ 2 \ 2 M 2 c 
M 2 -q 2 



(B.18) 



The dots refer, again, to contributions of excited states and different Lorentz structures that 
we do not consider. Expressing |Aj| 2 by the expression in Eq. (|B.10|) and Eq. (|E>. 1 1 
obtain the following sum rules: 



we 



n 



^M* B E 3 (s /Ml) + h -MlE x {s G /Ml) + £ U - ^ 

d \ B 



£ Li 
J I 



M B E 3 (s /M B ) + b -Ml El (s /Ml) + y (4 - ^| 

±M B E 3 (s /Ml) + b -M BEl (s /Ml) + 
M e B E 3 (s /M B ) + -MlE^so/Ml) + - U - ^ 

\ B 



O At. 

-M B E 3 (s /M 2 B ) - -MlE^/Ml) 

ft = ~ 1 • (B.19) 

QM B E 3 (s /M 2 B ) + -MlExisv/Ml) 

Inserting the numerical values of parameters, we end up with the following estimates for the 
higher-twist matrix elements: 

ff(n = 1 GeV) = 0.6 ± 0.2 , 
/2V = 1 GeV) = 0.15 ±0.06, 

= 1 GeV) = 0.22 ± 0.15 . (B.20) 



C Handbook of nucleon distribution amplitudes 

In this Appendix we give a complete list of all 24 nucleon distribution amplitudes. 

The chiral structure of our basic set of eight independent distribution amplitudes is 
summarized in Table 2 in the text. The corresponding representations for the remaining 
16 distributions are collected in Table [A] and Table [B[ The set in Table |A] contains the 
amplitudes that can be directly obtained from Eq. ( |2.26| ) to Eq. (|2.29|) by flipping spin- 
projections of the two up-quarks. Table || contains four distributions that are obtainable 
from the previous ones by exchanging light-cone projections, and also Ti,T 2 ,T 5 and T 6 that 
cannot be read off from the previous expressions in a straightforward way and therefore are 
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Lorentz-structure 


Light-cone projection 


nomenclature 


twist-3 






[VI + Ax] (a;,) 


twist-4 






[V 2 + A 2 \ ( Xi ) 








[V 3 + A 3 ] ( Xi ) 








[T 3 + T 7 + S 1 - P x \ ( Xi ) 


twist-5 




J 1 


[V 5 + As] (x^ 








[V* + M] (xi) 




(C jff)®]f 




[T 4 + T 8 + S 2 -P 2 ] ( Xi ) 


twist-6 






[V 6 + A 6 ] ( Xi ) 



Table A: Eight nucleon distribution amplitudes that can be obtained from Tab. ([I]) by flipping 
the spin projections of the two up quarks. 

given below: 

(0| e^fulia^Cia^uliazz)^ ^ k {a 3 z) \P) = -2pz ^ Jv xe - ipz ^ x ^ Tfa) , 
(0\e ijk (ul(a 1 z)Cia ±z u l j (a2z) S j 7 X jd l k (a 3 z) \P) = -2pz ^ Jvx e~ ipz ^ x ^ T 2 { Xi ) , 
(0| e ijk (ui(a 1 z)Cia ±p u i j (a 2 z)) 7^ ?4(a 3 z) \ P ) = " M ' ^ Jvxe~ vpz ^ x ^ T 5 { Xi ) , 

(0| e ijk (ul(a lZ )Cia ±p u}(a 2 z)^ $d l k {a 3 z) \P) = —M 2 ^ Jvxe~ ipz ^ x ^ T 6 (x 4 ) .(C.l) 
Note also that for scalar and tensor contributions the following representations can be given: 
(0\e^ k (ul(a lZ )Cu}(a 2 z)^dl(a 3Z )\P) = ~M frf Jvx e -**S>*«* [S 1 + P 1 ] ( Xi ) , 

(0|^ fe (wJ(a 1 z)Cz ( 7^t(a 2 z))/4(a3z)|P) = ^pzM ^ Jv X e~ ipz ^ x ^ [T 3 -T 7 ] ( Xi ) . 

(C.2) 

In the following subsections we give the complete set of distribution amplitudes starting 
from twist-3 to twist-6. The numerical values of the expansion parameters can be obtained 
from Table ^| in the text. 

C.l Twist-3 distribution amplitudes 

Vi(xi,/i) = 120 Xl x 2 x 3 [</>l{fj) + 03 - 3x 3 )] , 
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Lorentz-structure 


Light-cone projection 


nomenclature 


twist-3 


(Cia ±z ) ® 7 V 


ti-j- li-j- 


Ti(xi) 


twist-4 


(Cia ±z ) <g> 7 X $ 




T 2 {xi) 








-[T 3 -T 7 -3 1 -P 1 ] ( Xi ) 








- [T 3 + T 7 - Si + P x \ fc) 


twist-5 


{Cia ±p ) <g> 7 X / 


u l u l d i 








u+ u+ d-^ 


-[T 4 -T 8 -S 2 -P 2 ] (a*) 








-[T 4 + T 8 -^ 2 + P 2 ] ( Xi ) 


twist-6 


(Cia ±p ) (g) 7 X ^ 




T & {xi) 



Table B: The remaining eight nucleon distribution amplitudes that enter the expansion in 
Eq. (p|). 



Ax(xi, fi) = UOx^x^x-i - x x )(j) 3 (/x) 



Ti(xi,/j,) = 120x!X 2 x 3 
C.2 Twist-4 distribution amplitudes 



(C.3) 



= 24xix 2 [0°O) + 0i" - 5x 3 )] , 
A 2 (xi,/j,) = 2Axix 2 (x 2 -xi)04 (/i) , 

VsO^aO = 12x 3 [^2(a*)(1 - x s) +^{^){x\ + x\ - z 3 (l -x 3 )) +^ 4 f (/i)(l -x 3 - 10xix 2 )] , 

A 3 (xi,ii) = I2x 3 {x 2 -x x ) [(ipl(fj) + +^(n){l — 2xs)] , 

T 3 ( Xi ,ti) = Qx 3 [(C + 0° + - x 3 ) + (C + <K ~ Tj>Z)(jj)(xl + x\- x 3 (l - x 3 )) 

+ (& + <Pt + V^)(a0(1 -xa- Wx 1 x 2 )] , 
r 7 (xi, n) = 6x 3 + 0° + ipl)(n)(l - x 3 ) + (-£ + 01 - ^DMfri + ^2 - zs(l - z 3 )) 

+ + 01 + ^)(^)(1 " ^3 - 10xiar 2 )] , 
T 2 (x,,/i) = 2437X2 [£ 4 °(/i) -5x 3 )] , 

SiOi, /i) = 6x3(22 - xi) [(£ + 0^ + V 4 ° + £ 4 + + 04 + ^4 + )(aO + (C + 01 - ^4~)0")(1 " 2a*)] , 

/i) = 6X 3 (X 2 - Xx) [($ - 04* - V>4 + ^4 + - 01 - ^4 + )(^) + (C ~ 01 + ^")(a0(1 " 2x 3 )] ■ 

(C.4) 
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C.3 Twist-5 distribution amplitudes 



V 5 { Xi ,n) = 6x 3 + (pt(fi)(l - 2x 3 )] , 

A 5 (xi,fx) = 6x 3 (x 2 - xi)05 (//) , 

^(xi, /i) = 3 [^s(^)(l - a*) + V> 5 ~(aO (2zix 2 - x 3 (l - x 3 )) + (A*)(l - a* - 2(x? + ^))] , 
A 4 (xi,/i) = 3(x 2 - xi) [-^5 (^) +i/j5(fj)x 3 +ip^(fi)(l - 2x 3 )] , 

Tfa, n) = I + ^5 + 0s) M(l - ^3) + fe" + 5 " - ^5") (^) (2^i^2 - ^(1 - x 3 )) 

+ (# + 5 + + ^ 5 + ) (A*)(l - a* - 2(x? + **))] , 
T 8 (x l? /.) = I + 0° - H) (A*)(l - x 3 ) + (05 - " £") (A*) (2*1*2 - x 3 (l - x 3 )) 

+ (0+ + ^ - e 5 + ) (A*)(l - x s - 2(x? + **))] , 
T 5 (x,, /i) = 6x3 $(//) + e 5 + (/^)(l " 2x 3 )] , 

S 2 (x,, /i) = |(X 2 - XO [- + 0° + £5°) (/i) + fe" + 05 - O 

+ (£+ + 0++^) M(1-2X S )] , 



P 2 (x,, /i) = -(X 2 - Xl ) M + ^ ~ e 5 °) (/X) + fe- - 05" + 



(/i)(l-2x 3 )] , 



(C.5) 



C.4 Twist-6 distribution amplitudes 



V 6 ( Xi , /i) = 2 [<f>°M + 6 + (/i)(l - 3x 3 )] 
Ae(xi,n) = 2(x 2 - xi)0e , 



T 6 (xi, fj) = 2 



4>°M + \(ti - 6 + )(l-3x 3 ) 



(C.6) 
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